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CN \ As the nearest dramatic example of a poststarburst galaxy driving a galactic 

\ wind, NGC 1569 is an ideal test environment to understand the impact of "feed- 

back" from massive star lives and deaths on the surrounding ISM. We present 

Q\ [ HST WFPC2 1 narrowband imagery of NGC 1569 in an attempt to understand 

the underlying ionizing emission mechanisms on a 3 pc scale and to generate a 
H II region size distribution and luminosity function. We use [0 lll]/H/5 and 
[S Il]/Ha ratio maps to find that non-photoionizing mechanisms (e.g. shocks) are 
responsible for 10±3% of the Ha emission, 2.3-3.3 times larger than results from 
Q-f similar galaxies. Note that our method of determining this result is different than 

these past results, a point that we discuss further in the article. The area of the 
non-photoionized region is 10-23% of the total. Our results for NGC 1569 indi- 
cate that these non-photoionized areas do not lie in low surface brightness regions 

.^J \ exclusively. A comparison with multiwavelength point source catalogs of NGC 

1569 indicates that the dominant non-photoionizing mechanisms are shocks from 
supernovae or winds from massive stars. To explain this large percentage of non- 
photoionized emission, we suggest that NGC 1569 is, indeed, in a post-starburst 
phase as previous authors have claimed. We also derived slopes for the H II 
region luminosity function (-1.00±0.08) and size distribution (-3.02±0.27). The 
luminosity slope, though shallow, is similar to previous work on this galaxy and 
other irregular galaxies. The size distribution slope is shallower than previous 
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slopes found for irregular galaxies, but our slope value fits into their confidence 
intervals and vice versa. Within 4 pc of the 10-20 Myr old superstar clusters Al, 
A2, and B, no bright HII regions exist to a luminosity limit of 2.95xl0 36 erg s -1 , 
suggesting that the winds and shocks have effectively terminated star formation 
in this small cavity. In the three annular regions around the SSCs, both the HII 
region luminosity function and HII region size distribution are consistent with 
respect to one another and the galaxy as a whole. The HII region surface den- 
sities within the annuli remain the same as the annuli are moved away from the 
superstar clusters. These results indicate that "feedback" effects in NGC 1569 
are confined to the immediate vicinity of the most recent massive star formation 
event on scales of ~ 1 pc. 

Subject headings: galaxies: individual (NGC 1569) — galaxies: ISM — galaxies: 
structure — ISM: structure 

1. INTRODUCTION 

1.1. Feedback Mechanisms 

Elevated star formation rates of starburst galaxies produce a corresponding increase 
in the numbers of massive stars and supernovae. Stellar winds and supernovae can ionize 
their surrounding interstellar medium (ISM) to levels that cannot be explained by radiative 
processes (photoionization) alone. These non-photoionizing (shorthand for shocks, etc.) 
mechanisms, therefore, may contribute an appreciable amount to the excitation of the warm- 
phase ISM. In recent studies of non-photoionizing mechanisms in starburst galaxies, Calzetti 
et al. (2004) found that 3-4% of the Ha emission in nearby starburst galaxies were produced 
by non-photoionizing sources and that the main cause of this non-photoionized emission was 
due to the mechanical energy injected from the recent and past star formation. These 
conclusions are significant because they are one of the first measurements of the small-scale 
feedback mechanisms from stars (i.e., stellar winds & supernovae explosions) produced in 
the starburst. Knowing the importance of star formation of the surrounding ISM on small- 
scales allows us to understand how larger amounts of star formation affect the global ISM, 
and thereby, galaxy evolution. Such global feedback mechanisms can drive the remaining gas, 
necessary to produce stars, from the galaxy to enrich the intergalactic medium (IGM). Such 
expansion of the ionized gas into the IGM appears to be occurring in NGC 1569, according 
to Ha (Heckman et al. 1995) and X-ray (Martin et al. 2002) studies. Thus, NGC 1569 is 
ideal for studying the current non-photoionized processes within the galaxy and determining 
how important these processes are to the future of the ISM composition in this galaxy. 
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1.2. H II Region Luminosity Function and Size Distribution 

The luminosity function (LF) of H II regions in a galaxy provides an excellent tracer of 
the recent (age < 10 Myr) star formation. The form of this function is typically given as a 
power law in the form: 

N(L)dL = AL~ a dL, (1) 

where N(L)dL is the number of H II regions with a Ho luminosity between L and L + dL, 
A is a constant, and a is the slope of the power law. Past studies (e.g., Kennicutt et al. 
1989) have found that the power law decreases for later galactic types (1.75 for LMC, 2.1 for 
M101, and 2.3 for M31). Some galaxies show a break in the LF slope where the fainter H II 
regions have a shallower slope (e.g., Rozas et al. 1996). Some suggest that this slope break 
is due to a transition between normal H II regions and giant H II regions (e.g., Kennicutt 
et al. 1989). Oey & Clarke (1998) suggest that the evolution of the ionizing stars in the 
clusters causes this slope break. 

The Ho luminosity function of NGC 1569 has been determined in the past by Young- 
blood & Hunter (1999). They find 31 H II regions and derive a slope of -1.38±0.16. They 
also find that the slope of the luminosity function turns over at lower luminosities. With the 
improved spatial resolution of HST, we find a factor of 30 more H II regions in NGC 1569 
than was previously reported and a slope similar to that of Youngblood & Hunter (1999). 

The size distribution of H II regions provides insight into the evolutionary state of star 
forming regions within a galaxy. This size distribution has been fitted by an exponential 
law in the past (e.g., van den Bergh 1981). However, Oey et al. (2003) suggest that the 
size distribution is best explained by a power law distribution because the LF and the 
size distribution are connected. For NGC 1569, Youngblood & Hunter (1999) fit the size 
distribution of the H II regions with an exponential function and derive a characteristic H II 
diameter of 34±3 pc. 

1.3. NGC 1569 

NGC 1569 is a nearby (D = 2.2±0.6 Mpc; Israel 1988) Im galaxy that has been well- 
studied over the last 20 years. Two prominent, stellar-like features are located at the center 
of this galaxy, which are thought to be super-star clusters (SSCs; see Prada et al. 1994) 
formed in a starburst 2-10 Myr ago (Gonzalez-Delgado et al. 1997). However, SSC A has 
been shown to be two stellar clusters superimposed (Maoz et al. 2001). Several studies 
using optical interference-filter imagery of the ionized gas show evidence of an eruptive event 
which occurred in the galaxy's past (e.g., Devost, Roy, & Drissen 1997). Kinematics of the 
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ionized gas that were studied by Tomita et al. (1994) showed that the expanding gas moves 
at speeds from 10 to 100 km s -1 . Heckman et al. (1995) found that the optical filaments at 
distances of 2 kpc from the center of the galaxy are traveling over 200 km s -1 . 

The consensus of all these and other past studies is that a starburst occurred approxi- 
mately 10 Myr ago. This event (of unknown origin but probably produced from an interaction 
with a H I companion, Miihle et al. 2005) produced super star clusters (SSCs) Al, A2, and 
B. Because of the numbers of massive stars and their rapid evolution (first giving rise to 
stellar winds and then supernovae), the galaxy underwent a pronounced kinematical and 
morphological change, even disruption, during the past several Myr. Evidence suggests that 
some fraction of the supernova ejected material will escape the galaxy and will enrich the 
IGM (Martin et al. 2002). 

With the superior angular resolution of the Hubble Space Telescope over ground-based 
instruments, we determine the current amount of non-photoionized Ha emission within the 
galaxy on a 3 pc by 3 pc scale. Using previous results, we determine the mechanisms 
necessary to produce this non-photoionized Ha emission in this galaxy and whether the 
past star formation can explain the current amount of non-photoionized emission within the 
galaxy. We have also generated a new H II region luminosity function and size distribution. 
The observations and data reduction are presented in §2. Basic results of our analysis 
specifying the fraction of non-photoionized gas in NGC 1569 are given in §3, and discussion 
of our findings is given in §4. Results and Discussion of the H II region size distribution and 
luminosity function are given in §5 and §6 respectively. In §7 a summary of our findings and 
concluding remarks are made. 



2. OBSERVATIONS AND DATA REDUCTIONS 

WFPC2 images of NGC 1569 were taken 1999 September 23 for the Cycle 8 program, 
GO-8133. NGC 1569 was oriented in two of the wide-field chips (WF2 & WF3) of the 
camera with a nearby 10"* magnitude star placed out of the field of view (see Figure 1). The 
effective plate scale is 0'.'0996 pixel -1 (1.06 pc pixel -1 at the adopted distance of 2.2 Mpc). 
The GO-8133 images used here were F486N (H/3), F502N ([O III]), F547M (~ Stromgren y), 
F656N (Ha), and F673N ([S n]). NGC 1569's heliocentric radial velocity is -104 km s -1 . 
The shift of each emission line is ~ —2 A, and therefore, all emission lines were observed very 
near the center of the filter transmission curve. Observational parameters for these data are 
found in Table 1. See Buckalew et al. (2000) for more details on the calibration processes 
used for these images. In that paper, the reduction process was detailed for the F502N image 
only. However, the reduction process is similar for the F486N (H/3), F502N ([O III]), F656N 
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(Ha), and F673N ([S n]) images. For the analysis of the non-photoionized gas fractions, the 
images were binned in 3 pixels x 3 pixels (3x3), 5 pixels x 5 pixels (5x5), 7 pixels x 7 pixels 
(7x7), and 9 pixels x 9 pixels (9x9). Use of these binned images reduced the uncertainties 
in the emission-line ratio maps. For the analysis of the H II region luminosity function and 
size distribution, a calibrated, unbinned Ha image was used. 

The filter transmission curve of F656N covers Ha as well as [N n] AA6548,6583. To 
remove this contamination of the Ha emission, we compared the eight values of [N ll]/Ha 
and [N n]/[S n] from Buckalew et al. (2005) to determine which varied the least. Calzetti et 
al. (2004) determined that [N n]/[S n] varies little over similar starburst galaxies. However, 
we found that [N ll]/Ha varies less in NGC 1569 than [N n]/[S n]. Taking into account 
the blueshift, the wavelengths have changed by 2 A for the 3 emission lines. The total 
system throughput ratios of A6548/A6563 and A6583/A6563 were determined, and the total 
percentage of [N n] contamination is found to be 1%. The Ha image was multiplied by 0.99 
to reflect this amount of contamination. 

Flux calibration of the continuum-subtracted, interference-filter images was done using 
the procedure from Buckalew et al. (2000). The flux conversion numbers are in units of erg 
cm" 2 DN" 1 and are 1.35 x 10~ 14 for F487N, 1.07 x 10^ 14 for F502N, 4.14 x 10~ 15 for F656N, 
and 3.95 x 10~ 15 for F673N. 

The image with the lowest signal-to-noise, H/3, was used to generate a mask. We deter- 
mined the standard deviation of the background in the H/3 image at several sky locations. 
The standard deviation was 2.5xl0 -17 erg s _1 cm -2 in the 3x3 H/3 image. Standard devia- 
tions were also measured for the 5x5, 7x7, and 9x9 binned images. We generated H(3 masks 
at the 3a, 5a, 7a, and 9a levels above the mean sky to multiply to the appropriate binned 
images. The results discussed below will focus predominantly on the 3x3 binned images 
that have been masked to remove pixels with a flux level less than 5a above the mean sky 
(hereafter 3x3 5a). 

Three ratio maps were made of the galaxy: Ha/H/3, [O lll]/H/3, and [S ll]/Ha. The 
Ha/H/3 ratio map was used to deredden the [O lll]/H/3 and [S ll]/Ha maps. The pixels 
in the Ha/H/5 ratio map with Ha/H/3 < 2.86 were removed because 2.86 is our assumed 
theoretical ratio of these two hydrogen lines (Osterbrock 1989), typical for a H II region 
with T c of 10,000 K and N e of 100 e~ cm -3 . The following equation was used to convert 
Ha/H/3 to C(H/3): 

C{Hf3) = (log[Ha/Hf3] - Zo#[2.86])/0.326, (2) 
where Ha/H/5 indicates a value from our ratio map. With the C(H/3) ratio map, we dered- 
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dened the other ratio maps using the following formulae: 



I([0 III])/I(Hf3) = F([0 III])/F(Hf3) x i(T 003 W) ? 
I([S II])/I{Ha) = F([S II])/F{Ha) x nr - 020 ^, 



(3) 



(4) 



where IQ means the dereddened fluxes and F() indicates the reddened fluxes from the ratio 
maps. The multiplicative factors to C(H/3) are based on the reddening law from Seaton 
(1979). The average C(H/3) for this galaxy is 0.69 (E(B-V) = 0.47 & A v = 1.5 mag). 
Thus, the reddening only decreases the ratios by an average of 5%. With the dereddened 
ratio maps, we determined which 3x3 5a [O in]/H/5 and [S n]/Ha ratios had a reasonable 
Ha/H/3 ratio (i.e. Ha/H/3 > 2.86). Figure 2 displays the 3x3 5a [O ni]/H/3 ratio map, and 
Figure 3 displays the 3x3 5a [S irj/Hct ratio map. 

3. EMISSION MECHANISM RESULTS 

3.1. Determination of Photoionized and Non-Photoionized Regions 

In Figure 4 we plot the 3x3 5a [O lll]/H/3 versus [S irj/Hct values to compare with 
the maximum starburst relation from Kewley et al. (2001). This line denotes the upper 
limit between photoionized gas (below the line in Figure 4) and ionized gas from shocks or 
other non-photoionizing sources (above the line in Figure 4) . All points above this maximum 
starburst line are counted as emission due to non-photoionizing processes. All points below 
the maximum starburst line are counted as due to photoionizing processes. The distribution 
of the NGC 1569 points is similar to that of NGC 4214 (Calzetti et al. 2004). From Figure 
4 we find that 16% of the total number of points fall above the maximum starburst line. 
This areal percentage is similar to those found for NGC 3077, NGC 4214, NGC 5236, and 
NGC 5253 (Calzetti et al. 2004). The ratio of the Ha luminosity in these non-photoionized 
points to the total 3x3 5a Ha luminosity is 15%. 

Compared to previous results from Calzetti et al. (2004), this percentage is four to 
five times larger than the percentages of NGC 3077, NGC 4214, NGC 5236, and NGC 5253 
(range 3-4%). Note that the method for determining non-photoionized results is different 
between ours and that of Calzetti et al. (2004). We choose everything above the maximum 
starburst line while Calzetti et al. (2004) use an additional criterion (their dashed line in 
their Figure 7). This dahsed line, in a basic way, limits log([Sn]/Ha) to values >-0.5 or >- 
0.4. These differences in selection criterion are based on differences in the shock models used. 
We use the shock models of Hartigan et al. (1987); see Calzetti et al. (2004) for their shock 
models. If we use their additional limitation, we determine a non-photoionized percentage of 
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2-1%. Thus, we can simply say that the non-photoionized percentage of emission is different 
in NGC 1569 than that of other irregular galaxies regardless of the shock models used. 

The 3x3 3a, 3x3 7a, and 3x3 9a non-photoionized Ha luminosity percentages of 
NGC 1569 are 19%, 12%, and 10% respectively. The limited decrease in non-photoionized 
percentage as a function of increasing a suggests that the non-photoionized emission of 
NGC 1569 is found in high surface-brightness regions. This lack of decrease is at odds 
with what Calzetti et al. (2004) found for other, similar irregular galaxies. The fraction 
of non-photoionized flux found in the low signal-to-noise threshold (5<r) maps should be 
much larger (not just a factor of 1.5) than the fraction found in the higher signal-to-noise 
threshold maps. The lower signal-to-noise ratio maps should give a higher percentage simply 
because more lower surface brightness regions would be counted, typically the location of 
non-photoionized emission. We discuss in §4 possible sources to explain this high surface 
brightness non-photoionized emission in NGC 1569. 

Because the percentage found from the 3x3 5a maps is so much larger than the previous 
results, we determined the fraction of non-photoionized Ha emission from the 5x5, 7x7, and 
9x9 5a maps. The 5a non-photoionized Ha luminosity percentages for 5x5, 7x7, and 9x9 
are 16%, 10%, and 10% respectively. Increasing either the bin size to 9x9 or the sigma to 
9a causes a 5% decrease in the non-photoionized Ha luminosity. The 9x9 9a result is 5%, a 
factor of 3 decrease, but still the largest percentage for an irregular galaxy. The majority of 
these results fall at 10%±3%, which is the percentage that we claim for NGC 1569, whether 
we use 3x3 5a or 9x9 9a maps. The decrease in the fraction of non-photoionized luminosity 
with bin size is caused by spatial dilution with the photoionized emission. 

The idea of spatial dilution brings up the important point that NGC 1569 is closer, by 
a factor of 2, than the irregular galaxies studied in Calzetti et al. (2004). Since they use a 
3x3 5a result, we must compare a similar spatial result to theirs. The results above that are 
comparable in spatial resolution are the 5x5 5a and 7x7 5a results (the best being a 6x6 
5a result). Both of these results are 10%. Thus, our final, accepted percentage of 10%±3% 
works well for comparison with these further irregular galaxies. 

3.2. Further Quantities Derived from the Ha Images and Ratio Maps 

From the Ha luminosity of the 3x3 5a results, we derive the star formation rate (SFR) 
and number of hydrogen ionizing photons per second (Q ) using the equations from Kennicutt 
(1998). These results, along with those from §3.1, are summarized in Table 2. When we 
compare these results to those of similar nearby starburst galaxies in Calzetti et al. (2004), 
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we find that the SFR and Q of NGC 1569 are similar. This similarity is also true of the 
total Ha luminosities and areal coverage of the non-photoionizing luminosity. However, we 
showed in §3.1 that the major difference is the ratio of non-photoionized Ha luminosity to 
the total luminosity, which we will discuss further in §4. 

We have compared the properties of the photoionized and non-photoionized properties 
using a Wilcoxon test to determine differences between the samples. As confirmation that 
the Wilcoxon test works, we find that [O lll]/H/3 and [S ll]/Ha are much higher in non- 
photoionized areas compared to photoionized areas. We find that the southern areas of NGC 
1569 are prevalently non-photoionized areas. The visible outflow and Ha shell lie on this side 
of NGC 1569. This outflow lies on the near side of NGC 1569, making the non-photoionized 
emission more prominent. We also find that non-photoionized areas tend to have larger 
C(H,3), which can be explained by dust or non-photoionized mechanisms. Hartigan et al. 
(1987) shows that shock excitation can elevate Ha/H/3. Finally, we find that the continuum 
flux is higher at the locations of the non-photoionized emission compared to the locations of 
the photoionized emission. Perhaps, this finding means that the non-photoionized emission 
is typically coming from areas closer to stars and star clusters, but this interpretation is 
speculative at best. What we find below is that while some non-photoionized emission 
is near bright continuum sources, the majority is scattered throughout the galaxy neither 
tending to be near or far from resolved stars and star clusters. 



4. DISCUSSIOIN OF NON-PHOTOIONIZED EMISSION 

A significant fraction of the Ha luminosity in this galaxy arises from non-photoionizing 
sources. NGC 1569 is riddled with Ha shells, arcs, and filaments throughout the main 
body of the galaxy (see Figure 1). Can the mechanical energy from the stellar winds and 
supernovae explain the amount of non-photoionized Ha luminosity? Do potential sources 
such as Wolf-Rayet stars and supernovae remnants or different gas phases, such as the X-ray 
hot gas, explain the locations of some of these non-photoionized sources? 



4.1. Is the Starburst Enough? 

4-1.1. Continuous Star Formation 

We derive a SFR of 0.15 M yr -1 from the Ha luminosity. We generated Starburst 
99 (Leitherer et al. 1999) continuous star formation models from Myr to 100 Myr to 
determine the amount of mechanical luminosity generated over this time frame. The model 
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parameters are the metallicity of Z = 0.004, lower and upper stellar mass limits of 0.1 M 
and 100 M , and a Salpeter initial mass function (IMF). The 10 Myr and 100 Myr model 
mechanical luminosities are 3.02 xlO 40 ergs s _1 and 8.91 xlO 40 ergs s _1 respectively. We 
follow the prescription of Calzetti et al. (2004) and assume that the Ha luminosity due 
to the mechanical luminosity is 2.5%. The Ha luminosity due to these mechanical energy 
sources is 7.5xl0 38 ergs s _1 and 2.2xl0 39 ergs s^ 1 for 10 Myr and 100 Myr respectively. 
Our estimate of the non-photoionized Ha luminosity is 2.9xl0 39 ergs s _1 . These model 
non-photoionized Ha luminosities are 4-1.3 times too low to explain the current amount of 
non-photoionized Ha luminosity that we estimate. Greggio et al. (1998) state that ~5 Myr 
ago the SFR was 0.5 M yr -1 and that this SFR was constant over the past 100 Myr. If we 
assume that the SFR was 0.15 M yr -1 over the past 5 Myr, the mechanical luminosity of 
this contribution is 10 times too low to explain the non-photoionized Ha luminosity. The 
SFR of 0.5 M yr~ x can provide a sufficient amount of mechanical luminosity after ~10-20 
Myr to explain the current non-photoionized Ha luminosity seen. We agree with Greggio et 
al. (1998) in so far that the SFR must have been higher in the past to explain the amount of 
non-photoionized Ha emission currently seen, if this emission is from the mechanical energy 
of supernovae and stellar winds. 



4-1.2. Instantaneous Star Formation 

We assume that the star formation occurred instantaneously and that the stellar mass 
produced in the burst is equivalent to the most massive star clusters. According to Anders et 
al. (2004), the mass of the three largest clusters is between 6.2-6.6 log(M ); similar masses 
were found in Buckalew et al. (2005). We used the Anders et al. (2004) mass range in 
an instantaneous burst Starburst 99 model (other parameters are similar to those used in 
the continuous star formation) to determine the lower and upper mechanical luminosities at 
10 Myr and 100 Myr, which were converted to Ha luminosities using the same assumptions 
as before. The Ha mechanical luminosity range for 10 Myr and 100 Myr is 7.1xl0 38-39 
ergs s _1 and 1.6xl0 32-33 ergs s -1 respectively. At 100 Myr, the non-photoionized Ha lu- 
minosity remaining long after the burst is minimal. But at 10 Myr, the non-photoionized 
Ha luminosity is similar to that found in our study, and a majority of the mass range gives 
enough non-photoionized Ha luminosity to explain our results as well. This result is con- 
sistent with independent measurements of the star formation history that indicates that the 
burst probably occurred ~10 Myr ago (e.g., Greggio et al. 1998). It also suggests that the 
SSCs and bright H II region are responsible for the majority, if not all, non-photoionized Ha 
emission within this galaxy. However, not all non-photoionized Ha luminosity is necessarily 
explained by these sources. Some non-photoionized emission may be explained by individual 
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events/sources. 

4.2. Sources Responsible for the Morphology of the Non-Photoionized Gas 

We have already shown that the past star formation can explain the amount of non- 
photoionized gas. However, we can explicitly show that certain types of objects (e.g. Wolf- 
Rayet stars) are responsible for specific non-photoionized emission and the morphology of 
that emission. We compare the locations of the non-photoionized Ha with the locations of 
detected radio supernovae, Wolf-Rayet stars, star clusters, radio point sources, X-ray point 
sources, X-ray extended emission, and CO molecular cloud detections. 

4-2.1. Wolf-Rayet Stars & Star Clusters with Wolf-Rayet Stars 

Wolf-Rayet stars make good sources of mechanical ionization because they have strong 
winds. They are also indicators that star formation has occurred recently. In Figure 5, we 
plot the positions of Wolf-Rayet stars from Buckalew et al. (2000) with the non-photoionized 
points on the Ha image. The positions of the Wolf-Rayet stars are 2 pixel-wide black squares 
and labeled with their name from Buckalew et al. (2000). The non-photoionized points are 
1 pixel-wide white squares. 

Wolf-Rayet star SI has a few non-photoionized points adjacent to its location. Wolf- 
Rayet stars S2 through S5 and Unknown He II sources Ul through U3 have very little 
non-photoionized gas surrounding them. An arc lies above and to the left of S5 in Figure 
5, but we would not expect this particular arc to be caused by S5. SSC A (C4) is probably 
the main cause of this arc. Wolf-Rayet star S6 shows 3 non-photoionized regions adjacent to 
or on top of its location. Wolf-Rayet star S7 has no associated non-photoionized emission. 
This non-detection may result from low signal-to-noise in the H/5 image, or the possibility 
exists that the wind from S7 is too slow to produce significant mechanical luminosity. 

In the star clusters with Wolf-Rayet stars, cluster CI has a few non-photoionized points 
adjacent to it. Buckalew et al. (2000) showed that this "Wolf-Rayet" star cluster also 
had extended He II emission, which is coincident with the non-photoionized points. This 
extended He II emission signifies a shock velocity of ~120 km s -1 , the shock velocity where a 
maximum of He II emission is produced (Garnett et al. 1991). The wall of non-photoionized 
points between CI and C3 could be due to a wind-wind interaction between these clusters. 
C2 has one point of shocked gas. This non-photoionized point lies in the same direction as 
the extended He II emission found in Buckalew et al. (2000), implying a shock velocity of 
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~120 km s" 1 . C3 has a significant amount of non-photoionized gas on top of and around this 
star cluster. The line of points running through it do not trace any noticeable arc or filament 
in the Ha image (Figure 5). C4 (aka SSC A) is adjacent to S5 and the non-photoionized 
arc near C4 is probably caused by stellar winds from C4. C5 has an arc of non-photoionized 
emission to its lower right. This young cluster is probably the culprit for the majority of 
this non-photoionized emission. However, the massive star clusters in that H II region may 
also be contributing significantly, possibly this is a wind-wind interaction between the two 
star clusters. The average distance between the line and cluster is ~6 pc. If we assume a 
shock velocity of 100 km s" 1 , then the shocked gas originated a minimum of 10 4 ~ 5 yrs ago. 
Given the brevity of the Wolf-Rayet phase, this time scale would suggest that the onset of 
the Wolf-Rayet phase is responsible for the line of non-photoionized emission. 



4.2.2. Star Clusters 

We have plotted the WFPC2 F555W images of GO-6111 (Greggio et al. 1998) and GO- 
6423 (Hunter et al. 2000) along with the non-photoionized points in Figure 6. The positions 
of important star clusters are designated by their name from Hunter et al. (2000) and a 
line connects that name to the approrpiate star cluster. The non-photoionized points are 1 
pixel-wide white squares. Four arcs of non-photoionized gas near SSC A and SSC B are very 
prominent. The arc closest to SSC A appears to be wrapped around a cluster/star south of 
SSC A. Perhaps the wind of this star is interacting with a nearby, denser ISM. However, this 
arc is apparent in the left panel of Figure 6 but not the right panel. This arc (although more 
a line) is apparent in Figure 1 and must not, therefore, be an artifact due to the coordinate 
transformation. The cluster/star is the nearest stellar object inside this arc. No other point 
source or extended source can explain this arc. The other arc is just west of SSC B. This 
arc is apparent in Figures 1 and 6 and could be due to a wind-wind interaction between SSC 
B and cluster 28 (designation from Hunter et al. 2000) and the other cluster/star south of 
28 that has no designation. The last arc is west of SSC A and is associated with cluster 10 
and possibly clusters 13, 14, & 15. The arc is between cluster 10 (a.k.a. CI from §4.2.1) and 
clusters 13-15, again a possible wind-wind interaction between these clusters. Another arc is 
present and associated with clusters 39 and 40. Cluster 39 is C5 from §4.2.1 and has already 
been discussed. Clusters 30 and 35 have non-photoionized emission associated with their 
position (see Figure 6). Clusters 29 and 34 have at least one non-photoionized point near 
their location and are bright X-ray sources that are mentioned in §4.2.4. Finally, a large shell 
of non-photoionized gas surrounds the clusters 6, 7, 9, and 10. This shell is supported by 
the winds of these clusters and extended X-ray emission (discussed in §4.2.4). Several other 
sources such as the thermal radio source M-l from Greve et al. (2002) and the molecular 
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cloud 3 from Taylor et al. (1999) are also present at this region and help to show that this 
location is currently a site of vigorous star formation. 

4-2.3. Radio Sources from Greve et al. (2002) 

In Figure 7 we plot the radio-detected sources, such as supernovae and supernova rem- 
nants (SNRs), from Greve et al. (2002) along with the non-photoionized locations. The 
positions of the radio sources are 2 pixel-wide black squares and labeled with their name 
from Greve et al. (2002). The non-photoionized points are 1 pixel-wide white squares. The 
grayscale scheme used in §4.2.1 is implemented here. Note that radio sources M-a through 
M-d are tentative detections and are not known to be radio supernovae or SNRs yet. 

M-l and M-b,c,d are found inside a large shell of non-photoionized points. These points 
surround the brightest H II region complex within the galaxy. M-l, a thermal source ac- 
cording to Greve et al. (2002), is probably just the thermal radio emission from this large 
H II complex termed region 2 in Waller (1991). M-2 is a radio supernovae or SNR with 
one non-photoionized point adjacent to our marker. M-3 is a radio supernovae or SNR ac- 
cording to Greve et al. (2002). This source is situated in the middle of a large group of 
non-photoionized points. M-3 is more than likely responsible for some of these. However, 
we show in the next section that a X-ray binary is also found in the same area. According 
to Greve et al. (2002), M-4 is a thermal source (i.e., the H II region found below the point 
marker), and a few non-photoionized points lie near the marker. M-5 is a thermal source, 
according to Greve et al. (2002) and is probably associated with the coincident H II region. 
However, a Ha shell is found at the upper right hand corner of this H II region and is the 
source of the large number of non-photoionized points to the upper right of M-5. We think 
that M-5 is this shell and not the H II region. M-6 is a SNR according to Greve et al. (2002) 
and was also detected optically by Shopbell et al. (2000). As a side note, the size of the 
SNR in the radio is smaller than the optical counterpart (17 pc versus 21 pc). The brightest 
section of the SNR is coincident with the non-photoionized points. VLA-8 is most likely a 
SNR as well as VLA-16 (Greve et al. 2002). Regardless of their definition, these objects 
appear to contribute to the non-photoionized emission in their surroundings. The only SNR 
or radio supernovae without any nearby non-photoionized gas is VLA-11, probably caused 
by the lack of signal-to- noise in our Ha image at the position of VLA-11. 
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4-2.4- X-ray Point Sources 

In Figure 8 we plot the X-ray point sources from Martin et al. (2002) along with the 
non-photoionized locations on the Ha image. The positions of the X-ray sources are 2 pixel- 
wide black squares. The non-photoionized points are 1 pixel-wide white squares. We label 
each Martin et al. (2002) source with their number from that article along with a one letter 
identification indicating the nature of the source: C for cluster, S for SNR, and X for X-ray 
binary. Objects 19C and 22C are clusters 29 and 34 from Hunter et al. (2000). Both clusters 
have at least one non-photoionized point adjacent to their position. More than likely, these 
points are associated with those clusters, especially 19C since no other object exists in this 
evacuated cavity. Martin et al. (2002) were unsure if 28S was indeed a SNR. Comparing 
the Greve et al. (2002) results to these, we find that 28S is M-6, which is also the optical 
SNR discovered by Shopbell et al. (2000). We discussed the non-photoionized gas around 
this point in §4.2.3. All X-ray binaries except for 21X, 25X, 26X, and 29X have some nearby 
non-photoionized gas. Of these without non-photoionized gas, 21X and 26X sit inside large 
shells which have shocked gas along their edges. Perhaps these two X-ray binaries or their 
predecessors are responsible for these shells. Two of the remaining X-ray binaries deserve 
comment. 14X lies near M-3, the SNR. Thus, the SNR may not be completely responsible 
for this large non-photoionized area, or 14X was misdiagnosed as a X-ray binary and is 
actually the X-ray detection of M-3. Similarly, 24X sits atop the line previously explained 
by winds from C5, a cluster with detected Wolf-Rayet stars, and is possibly responsible for 
some percentage of the mechanical Ha luminosity found at this location. 

4-2.5. Extended X-ray Emission 

In Figure 9, we show the 0.3-7 keV X-ray image from Martin et al. (2002) along with the 
locations of the non-photoionized gas. One correlation is seen between the extended emission 
and non-photoionized gas, the large X-ray extended emission located over the brightest H II 
region complex (marked with a black circle in Figure 9). This location was also detected as 
M-l in Greve et al. (2002), a thermal source attributed to the H II region complex termed 
2 in Waller (1991). 

4-2.6. CO Emission from Molecular Clouds 

We have plotted in Figure 10 the non-photoionized points on a CO map from Taylor et 
al. (1999). The different molecular clouds are labeled using the designations from Taylor 
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et al. (1999). The only clouds found near non-photoionizing sources are 3 and 4. The 
molecular cloud 3 is found in the extended X-ray and thermal radio emission associated 
with the brightest H II region complex. The cloud appears to be circumscribed by the non- 
photoionized gas. Also, the local minimum in CO emission between 3 and 1+2 is coincident 
with one side of the non-photoionized shell circumscribing 3. The molecular cloud 4 is found 
coincident with some non-photoionized gas but too few non-photoionized points are present 
in that location to define the morphology of the cloud. Other low level features not labeled 
on the CO map are noise. 

4-2.7. Summary 

The most important correlation of the non-photoionized points with individual sources 
is with the largest H II complex. A ring of non-photoionized emission encircles the H II 
region complex 2 Waller (1991). This complex is the largest site of current star formation in 
NGC 1569. This site is also host to a large thermal radio source, extended X-ray emission, 
molecular cloud cores, supernovae, Wolf-Rayet stars, and star clusters. The stellar winds and 
supernovae explosions caused by this vigorous star formation are causing a large evacuation 
of gas which is striking the nearby ISM. Future kinematical evidence may show that the 
velocity of the gas is greater than the escape velocity. If it will escape, this complex may 
be similar (yet smaller) to the H II complex that contained SSCs Al/2 & B. In a few Myr, 
NGC 1569 may have another large Ha arm similar to the one currently ejected due to the 
super star clusters' creation. 

Other features of the non-photoionized emission appear to be related to SSCs Al/2 & 
B. A few arcs of non-photoionized emission seem to be near these two clusters and could 
be attributed to stellar winds emanating from these clusters and interacting with the winds 
from other nearby star clusters. Several other star clusters with Wolf-Rayet stars, without 
Wolf-Rayet stars, with X-ray emission, and/or with thermal radio emission appear with non- 
photoionized emission on top of or adjacent to these systems. The most important is C5 
(Buckalew et al. 2000). A significant arc of non-photoionized emission appears around this 
star cluster with Wolf-Rayet stars. The onset of the Wolf-Rayet phase in the OB stars is the 
likeliest explanation for the shocked emission. 

Other point sources that have non-photoionized emission are radio and optical super- 
novae, low-mass X-ray binaries, and Wolf-Rayet stars. With some of Wolf-Rayet stars, we 
notice that the position of the non-photoionized emission is in the same direction as the 
extended He II emission emanating from these sources. We suggest that this indicates the 
shock velocity is ~120 km s _1 , the velocity where He II 4686 A emission peaks. 
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4.3. Origin of the large percentage of non-photoionized gas in NGC 1569 

Calzetti et al. (2004) stated that the maximum percentage of non-photoionized Ha 
luminosity should be no more than 10-20%. NGC 1569 lies within this range. These max- 
imum percentages occur when the number of supernovae or stellar winds shocking the ISM 
stays the same and the level of photoionized emission is decreasing (i.e., when the OB stars 
are exploding as supernovae). This scenario typically occurs for older ages of a starburst, 
6-20 Myr roughly. Previous papers (e.g., Gonzalez-Delgado et al. 1997) concluded that 
NGC 1569 is in a post starburst phase. The ages of the super star clusters or large star 
clusters in NGC 1569 have ages around 10-20 Myr. Thus, the major star formation event is 
over, but star formation still occurs in smaller areas of the galaxies. The shocks produced 
from the major star formation event, which produced SSCs Al, A2, and B, are producing 
the non-photoionized Ha emission while a new, smaller (by a factor of 3) amount of star 
formation provides the photoionized component. 



5. MEASURING THE H II SIZE DISTRIBUTION AND LUMINOSITY 

FUNCTIONS 

We employed Hllphot (Thilker et al. 2000) to determine the positions, sizes, and fluxes 
of H II regions in NGC 1569. We used the values of 0.75, 1, 1.25, 1.5, 2, 3, and 10 for the input 
terminal gradient parameter that determines when the code stops determining the boundary 
of the H II region. The program generated ~6500 potential H II regions over the entire galaxy. 
Because the galaxy has a tremendous filamentary structure, the program interpreted each 
arc as a separate H II region. To remove the detections of the filamentary structures as H II 
regions, we accepted only those sources with a signal-to-noise of 100 or greater. After visual 
inspection to certify that the list does not contain filamentary structures, the actual number 
of H II regions used in the analysis was reduced to 1018. To deredden the luminosities of each 
source, we took the positions of the H II regions and calculated the mean 3x3 5a Ha/H/3 
value at these positions. The aperture used to calculate this mean was equal to the size 
of the H II region. In some cases, no valid points were found around the H II region, and 
we accepted a C(H/3) value equal to the average photoionized C(H/3) of the 3x3 5a results 
(0.67). The final range of Ha luminosities of the H II regions was io 36 - 4 -10 38 - 91 ergs s" 1 . The 
effective FWHM measurements of this final luminosity sample of H II sources were converted 
to 3a H II region diameters, producing a range of 6.3-44.8 pc. 
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5.1. Computing the H II Region Size Distributions 

A diameter histogram of the ~1000 H II regions was generated. The size of the bins 
are dictated by the size of the sample and by the standard deviation of the sample. The 
equation relating these three values is given by the equation of Scott (1979): 

Binsize = 3.5an~ 1//3 , (5) 

where a is the standard deviation and n is the sample size. These bins stretched to include 
the minimum value at the far left of the first bin and were added until the entire sample 
was covered. Uncertainties for each bin were estimated by assuming that the uncertainties 
of the bin were Poissonian (i.e. y/n). Because we plot the log of the bin sizes in Figure 11, 
the uncertainties of each bin are (-v/"ln[10]) _1 . 

Using the frequency values and uncertainties, we assumed that the H II region size 
distribution followed a power law like Oey et al. (2003). Linear regression was used to 
determine the slope (i.e. power law index) with the logarithmic values of the frequency (i.e. 
number of points in a bin) and diameter. The first fit taking into account all bins was 
significant to 23%, which is not a statistically significant fit. Removing the aberrant second 
left value from the fit, we derive a marginally significant (10%) fit to the data. The equation 
is the following: 

log(N) = -3.02 ± 0.27 x log(D) + 5.16 ± 0.27, (6) 

where N is the frequency of a certain bin and D is the diameter value of the middle of 
the bin. Others (e.g., Youngblood & Hunter 1999) have fit the size distribution with an 
exponential function. We have also fit the size distribution with an exponential function. 
The significance of such a fit (15%) is lower than the power law fit. For completeness, the 
equation is 

jy _ e -0.23±0.02xD+7.54±0.23 ^ 

where the letters represent the same values as the previous equation. 



5.2. Computing the H II Region Luminosity Functions 

Using the logarithmic Ha luminosities of the H II regions, we generated a luminosity his- 
togram. We follow the procedure outlined in the previous section for generating histograms. 
This histogram, displayed in Figure 12, was fitted by a power law distribution. The 3 left- 
most points are found below the fourth. The fit of all points except the four leftmost is 
significant to 8%: 

log(N) = -1.00 ± 0.08 x log(L) + 39.25 ± 2.14, (8) 
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where N is the number of points in one bin and L is the average luminosity for the bin. For 
comparison, the slope of the luminosity function for NGC 1569 in Youngblood & Hunter 
(1999) is -1.38±0.16, which is not within the 3a confidence interval of our result. The 
luminosity range of Youngblood & Hunter (1999) is between Log(L)=38.0-39.9. Our fitted 
logarithmic luminosity range is 36.4-38.91. We have plotted the slope of Youngblood & 
Hunter (1999) in Figure 12 over the luminosity range of their data. Note that their slope 
does work for our points found in a similar luminosity range. The most significant fit is for 
the luminosity range 37.22-38.91 and is the following: 

log(N) = -1.22 ± 0.06 x log(L) + 47.504 ± 2.14. (9) 

This slope is found within the 3a confidence intervals of the slope from Youngblood & Hunter 
(1999) and vice versa. The best fit equations along with the correlation coefficients for both 
the size distribution and luminosity function are given in Table 3. 

Because of the turnover of the luminosity function at the logarithmic luminosity value 
~37.0, incompleteness corrections may be necessary. We used IRAF to add elliptical galaxies 
to the image for HIIPHOT to detect. Elliptical galaxies were used because they are more 
extended than point sources and because no random H II region maker exists. 

The luminosity range of elliptical galaxies ranged from logarithmic values of 36.4 to 37.2. 
Above the logarithmic value of 37.1, HIIPHOT detects 100% of the added "H il" regions. 
Below this value, HIIPHOT fails to detect 8% of the added H II regions (incompleteness 
corrected bins are displayed in Figure 12). Adding this correction factor to those luminos- 
ity bins produces little change in the slope value. The slope value for the incompleteness 
corrected bins produces a slope change of 0.01. This small change is within the confidence 
interval of the original slope. Thus, the contribution of the incompleteness corrections is 
insignificant. With this incompleteness method, nothing can be said about corrections to 
the size distribution. The added sources are not H II regions but elliptical galaxies with 
similar Ha luminosities. 

We also produced luminosity functions and size distributions of three subsamples of 
the H II regions (see Figure 13). These three subsamples were determined based on their 
position from the "center of the starburst." We place the center of the starburst at the 
location halfway between SSCs A(l and 2) and B. The location of all H II regions in our 
sample have a minimum distance of 4 pc (maximum of 630 pc), suggesting that the winds and 
shocks have effectively terminated star formation in this small cavity. The first annuli is 115 
pc wide with the second is 62.7 pc wide, indicating that massive amounts of star formation 
are evacuated from the super star cluster area but are occuring rapidly outside this area. 
The luminosity function and size distribution fits for these subsamples were all significant, 
but the slopes for the size distribution and luminosity function are the same regardless of 
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position from the super star clusters. The H II region surface density was computed for each 
annuli, and these densities within the annuli increase by a factor of 1.33 from the first to the 
second annulus. These results indicate that these feedback effects in NGC 1569 are confined 
to the immediate vicinity of the most recent massive star formation event on scales of ~ 1 
pc. 

6. DISCUSSION OF H II REGION SIZE DISTRIBUTION AND 

LUMINOSITY FUNCTION 

6.1. H II Region Size Distribution 

Our power law slope of -3.02 is smaller than slopes of similar irregular galaxies (-3.39 to 
-4.16; Oey et al. 2003). Oey et al. (2003) state that a power law slope would provide only a 
good description of the size distribution for H II regions greater than 130 pc in diameter and 
that a slope should be zero for diameters less than 130 pc. All H II regions in our sample have 
diameters less than 130 pc, but we do not find a slope of zero. We speculate why we do not 
find a zero slope as either due to the fact that we have not detected several large (~45 pc) 
H II regions or that the resolution differences between our and past images are significant. 
The first possibility can be ruled out because our HST Ha images have exceptionally high 
signal-to-noise and high spatial resolution. Such large systems would be detected easily. 
Thus, the differences in spatial resolution of HST make this assertion untrue that the power 
law slope should be flat for diameters less than 130 pc. 

Comparing our number of H II regions with those from Youngblood & Hunter (1999), 
we have 30 times more H II regions. Our largest H II region diameter is 3 times smaller (45 
pc versus 159 pc) than that observed by Youngblood & Hunter (1999), demonstrating the 
superior resolution of the HST images. 

6.2. H II Region Luminosity Function 

Our H II region luminosity function slope, -1.00, found in this study is similar to those 
from previous studies of irregular galaxies (-0.60 to -2.67 Youngblood & Hunter 1999). Even 
with the resolution of the HST Ha image, we get a similar slope. The HST image allows 
us to resolve H II region complexes that are listed in past ground-based observations (e.g., 
Waller 1991) into several smaller H II regions. These large H II region complexes were used 
by Youngblood & Hunter (1999). Their slope is steeper and covers a smaller luminosity 
range. Their slope fits our luminosity data over their luminosity range of 38.0-39.0 dex. We 
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derive a similar slope of -1.22 over a larger luminosity range (37.22-38.91 dex) if we only fit 
our histogram bins before they begin to turnover. Thus, the slope of Youngblood & Hunter 
(1999) is accurate to luminosities down to logarithmic Ha luminosities of 37.22 dex. 

We can test the statement of Oey & Clarke (1998) that a power law slope used on lumi- 
nosities below L sat (10 38 - 5 ergs s -1 ) should be flat according to their models. We measured 
only five H II regions above this saturation limit, and we do not measure a flat slope for our 
luminosity function. Therefore, a nearly unsaturated (the number of stars in a H II region 
does not sample the IMF well) sample of H II regions produces a luminosity function with a 
shallow, though statistically significant, slope. Oey & Clarke (1998) also state that the slope 
of the luminosity function should break in two locations, once at the average luminosity of 
the luminosity function and once at the saturated luminosity value of 10 38 5 ergs s" 1 . We fit 
a linear regression to the bins in Figure 12 that are below and above the average luminosity 
(10 37 16 for our results) as well as the saturation limit. The slope below the average luminos- 
ity is 1.08±0.44, the slope above the average luminosity and below the saturated luminosity 
is -1.24±0.07, and the slope above the saturation limit is -1.55±0.34. The significances of 
these fits are 55%, 97%, and 91% respectively. The slope below the average luminosity is 
different than the slope above the average luminosity but not significantly. The other two 
fits are statistically significant, but the slopes are within each other's 3a confidence intervals. 
Since Oey & Clarke (1998) discuss the slopes by themselves, we do find that the slope above 
and below the average luminosity are statistically significant. We do not find that the slope 
for luminosities greater than the average luminosity is different that that for slopes greater 
than the saturation luminosity. 



6.3. Consequences of a Large Amount of Non-Photoionized Emission 

Our finding implies a large amount of non-photoionized gas in a small galaxy, and leads 
us to suspect previous derivations of reddening and metallicities based on ground-based data. 
Possibly the metallicity of this galaxy is in fact larger than previously thought. However, 
[O m]/H/3 and [S ll]/Ha alone cannot give accurate metallicity measurements because the 
emission-line ratios are degenerate in metallicity and ionization parameter (Kewley et al. 
2001). Reddening could be lower than previously thought because the Ha emission relative 
to H/9 is increased by shocks and other non-photoionizing mechanisms. We compare the 
C(H,3) averages of the non-photoionized and photoionized gas using the 3x3 5<r maps. The 
average for the non-photoionized gas is 0.76 and that of the photoionized points is 0.67. 
We test the significance in the differences in the averages using the F*-test for means. We 
calculate a p-value of 0, implying a significant difference. When the average C(R(3) is taken 
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for the entire galaxy, the value is 0.69. Thus, the presence of non-photoionized gas increases 
the derived C(H/3) by only 0.02 (or Ay of 0.04 mag) over the pure photoionized average. 
Because reddening is not affected significantly by a large percentage of non-photoionized 
emission, we hypothesize that the metallicity will not be altered by more than 0.2 dex. 

To determine if the large non-photoionized emission has any noticeable effect on the 
luminosity function and size distribution, we used Hllphot on the HST Ha images of NGC 
4214 using a similar procedure to that described in §5.2 and §5.1. The main difference 
between NGC 4214 and NGC 1569 is the amount of non-photoionized emission. No other 
differences are significantly different between these two galaxies, such as metallicity, star 
formation rate, etc. These Ha images are the ones used by Calzetti et al. (2004). We find a 
statistically significant luminosity function slope of -1.30±0.07 over a range of luminosities 
from 36.90 dex to 38.64 dex (sample size of 438 objects). The slope is steeper than the one 
found for NGC 1569 over a luminosity range of 36.81 dex to 38.91 dex, -1.00±0.08. Neither 
slope is found in the 3a confidence interval of the other, but their confidence intervals do 
overlap. However, our other slope for the range from 37.22 dex to 38.91 dex, -1.22±0.06, is 
well within the 3a confidence intervals of the other. We find a statistically significant size 
distribution slope of -2.37±0.2 over a range of diameters from 7.28 pc to 124 pc. This slope 
is shallower than the -3.39±1.94 slope from Oey et al. (2003), but our slope for NGC 4214 
is found well inside their confidence intervals (but not vice versa). The slope is also found 
outside the confidence intervals of our slope for NGC 1569 and vice versa. The maximum 
diameters are also a factor of 3 different, 45 pc for NGC 1569 and 124 pc for NGC 4214. 
The difference in size range is not caused by a difference in distances (only a factor of 1.3). 
In this comparison of slopes and ranges to the non-photoionized emission, we can speculate 
that the difference in mechanical luminosity, the only significantly different property between 
these two galaxies, could cause an overpressurization in the ISM or strip the ISM around 
the H II regions in NGC 1569 compared to NGC 4214. Either scenario could shrink the 
detected H II region diameters in NGC 1569 compared to NGC 4214. The lack of difference 
in the luminosity function slopes suggests that the scenario of overpressurization in the ISM 
is the more likely choice. If the material around the H II regions was being stripped, then 
we would see a distinct difference in the luminosity range and luminosity function slopes in 
NGC 1569 (smaller) compared to NGC 4214 (larger). This overpressurization scenario can 
be checked by determining the density of these H II regions. The average density of the H II 
regions in NGC 1569 should be larger than the average density of NGC 4214. 
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7. SUMMARY 

We have taken the WFPC2 narrowband imagery of NGC 1569 and generated Ha/H/3, 
[0 in]/H/5, and [S ll]/Hct ratio maps. We have determined the areas of non-photoionized 
and photoionized emission on a 3 pc by 3 pc scale. The non-photoionized Ha luminosity 
of NGC 1569 is 10%±3% of the total Ha luminosity. This value is approximately 2.5 to 
3 times larger than percentages found in similar starburst galaxies (Calzetti et al. 2004). 
However, we did point out that our method of determining the percentages are different from 
these past results. If we use a method more similar to theirs, we find that our percentage is 
1.5-2 times smaller than these previous results. Thus, we can conclude the non-photoionized 
percentage is different for this object compared to others. 

One-half to two-thirds of the non-photoionized emission originates in high surface bright- 
ness areas. We can explain the amount of non-photoionized Ha luminosity using an instan- 
taneous burst which uses a mass at least equal to the masses of the largest three star clusters 
in NGC 1569. Or we can explain the amount of non-photoionized Ha luminosity using a 
continuous star formation scenario with a SFR of 0.15 M Q yr _1 for 5 Myr and then the 
SFR from Greggio et al. (1998) for at least 10-20 Myr. We conclude that the main non- 
photoionizing mechanism is shocks mainly produced by supernovae or massive stellar winds. 
Individual sources responsible for these shocks include supernova remnants, x-ray binaries, 
and star clusters. The most prominent non-photoionized area is the arc/shell surrounding 
the brightest H II region complex. This site of vigorous star formation will likely continue 
to produce stars for the next several Myr. 

We have also used the Ha WFPC2 image to generate H II region size distribution and 
luminosity function. We find that a power law best describes the size distribution. The slope 
derived for the H II region size distribution (a = —3.02 ±0.27) is flatter than slopes for other 
irregular galaxies, but fits well within the confidence intervals of the past results and vice 
versa. Our luminosity function slope agrees with the result of Youngblood & Hunter (1999) 
over a logarithmic Ha luminosity range of 37.22-38.91 dex. We also derive a flatter slope, 
-1.00±0.08, that is within the norm for irregular galaxies and that is over the logarithmic 
luminosity range of 36.81-38.91 dex. 

The mechanical luminosity, which drives the non-photoionized emission, does not ad- 
versely affect the estimates of reddening or metallicity from ground-based spectroscopy by 
more than 3%. However, the slope of the size distribution as well as the range of sizes for 
NGC 1569 is different than that for NGC 4214. This finding coupled with the similarity in 
the luminosity function slopes and ranges for these galaxies suggests that overpressurization 
of the ISM from non-photoionized emission could cause these problems. Future density mea- 
surements of these H II regions could show explicitly whether this speculation is true. Within 
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4 pc of the 10-20 Myr old superstar clusters Al, A2, and B, no bright HII regions exist to a 
luminosity limit of 2.95 xlO 36 erg s -1 , suggesting that the winds and shocks have effectively 
terminated star formation in this small cavity. In the three annular regions around the 
SSCs, both the HII region luminosity function and HII region size distribution are consistent 
with respect to one another and the galaxy as a whole. The HII region surface densities 
within the annuli remain constant as the annuli are moved away from the superstar clusters. 
These results indicate that these feedback effects in NGC 1569 are confined to the immediate 
vicinity of the most recent massive star formation event on scales of ~ 1 pc. 

With these results, we add further evidence that a single burst event ~10 Myr ago 
created the super star clusters Al, A2, and B as well as the other star clusters responsible 
for the dominant event producing the current non-photoionized emission and properties of 
H II regions. 
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Fig. 1. — Ha Images of NGC 1569. Both images have been binned in 3x3 pixel bins and are 
displayed together on the same scale to show the locations of the non-photoionized points 
(white dots in the lower image). The upper image displays only the Ha image of NGC 1569. 
The lower Ha image depicts the locations of the non-photoionized gas, exceeding the 5a 
threshold, as white dots. Super Star Clusters Al, A2, and B are labeled to the right of the 
objects. The grayscale bars indicate the levels of the image and are in the units of 10 14 ergs 

— 1 — 2 

s cm . 

Fig. 2. — [O lll]/H/3 ratio map of NGC 1569. This ratio map has the same labeling and 
notation as Figure 1. Each point represents the 5a 3x3 [O lll]/H/3 binned results used in 
this study. The grayscale bar to the right indicates what numerical values that the various 
shades of grey represent. 

Fig. 3. — [S ll]/Ha ratio map of NGC 1569. This ratio map has the same labeling and 
notation as Figure 1. Each point represents the 5a 3x3 binned results used in this study. 
The grayscale bar to the right indicates what numerical values that the various shades of 
grey represent. 

Fig. A.— Plot of log([0 m]/H/3) vs. log([S ll]/Ha) for NGC 1569. Each point represents a 
3x3 bin of pixels (3.2 pc on a side) with a H/3 value 5a above the mean zeroed background 
of the H/3 image. The solid line represents the maximum starburst line from Kewley et al. 
(2001). The position and distribution of the points is similar to the ratios of NGC 4214 from 
Calzetti et al. (2004). However, the percentage of non-photoionized Ha emission is greater 
than the percentage of NGC 4214, by a factor of 3. 

Fig. 5. — Wolf-Rayet star locations superimposed on the NGC 1569 Ha image. This image 
has been binned in 3x3 pixel bins. Superimposed as white squares are the >5a locations of 
the non-photoionized gas within NGC 1569. Superimposed in gray circles with black outlines 
are the locations of Wolf-Rayet stars or star clusters with the presence of Wolf-Rayet stars. 
Each circle is labeled with its designation from Buckalew et al. (2000). Orientation of this 
figure is the same as that of Figure 1. Note the arc directly below C5, indicating that the 
wind from this cluster is possibly driving into the nearby ISM and creating a localized region 
of non-photoionized Ha emission. 
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Fig. 6. — F555W images of NGC 1569 taken from the archival data of GO-6111 (left) and 
GO-6423 (right). Orientation of these images are given by the compasses in the respective 
images. Superimposed in white are the 5a 3x3 locations of the non-photoionized gas within 
NGC 1569. The numbers are the star cluster designations from Hunter et al. (2000) and 
are discussed in §4.2.2. 

Fig. 7. — Supernova remnants labeled on the NGC 1569 Ha image. Superimposed in white 
are the 5a locations of the non-photoionized gas within NGC 1569. Superimposed as gray 
circles with black outlines are the locations of radio supernovae and supernovae remnants 
from Greve et al. (2002). Each circle is labeled with its designation from Greve et al. (2002). 
Note that most of these point sources have some nearby or superimposed non-photoionized 
gas (see §4.2.3 for more discussion). 

Fig. 8. — X-ray point sources labeled on the NGC 1569 Ha image. Superimposed in white 
are the 5a locations of the non-photoionized gas within NGC 1569. Superimposed as gray 
circles with black outlines are the locations of X-ray point sources from Martin et al. (2002). 
Each circle is labeled with its designation from Martin et al. (2002). Note that most of these 
point sources have some nearby or superimposed non-photoionized gas (see §4.2.4 for more 
discussion) . 

Fig. 9. — Broadband X-ray image from Martin et al. (2002). Superimposed in white are the 
3x3 5a locations of the non-photoionized gas within NGC 1569. North is up, and east is to 
the left. We have circled in gray the only obvious coincidence of extended X-ray emission 
and non-photoionized gas. This region is the brightest H II region complex in the galaxy 
and the location of the strong thermal radio source M-l of Greve et al. (2002). No other 
associations are apparent. 

Fig. 10.— CO map of NGC 1569 from Taylor et al. (1999). Superimposed in black are 
the 3x3 5a locations of the non-photoionized gas within NGC 1569. We have labeled the 
CO detections as designated in Taylor et al. (1999). North is up, and east is to the left. 
The only CO emission objects near non-photoionized points are 3 and 4 (designations from 
Taylor et al. 1999). The position of object 3 is the same as the extended X-ray emission 
circumscribed by non-photoionized points and as M-l, the large thermal radio source from 
Greve et al. (2002). 
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Fig. 11. — Histogram of H II region diameters with the fitted size distribution overplotted. 
The solid line is the fit as described by the equation in Table 3. The points are the bins as 
derived using the equation of Scott (1979). The x-axis uncertainties represent the size of the 
individual bins. The y-axis uncertainties are la in length. 



Fig. 12. — Ha luminosity histogram with the fitted luminosity function overplotted. The 
solid line is the fit of all luminosity bins and is described by the equation in Table 3. The 
long dashed line is the fit of the luminosity bins covering the range from 37 to 38.91 and is 
described in Table 3. The points are the bins as derived using the equation of Scott (1979). 
The x-axis uncertainties represent the size of the individual bins. The y-axis uncertainties 
are la in length. The dashed line in the bottom right corner is the slope as derived by 
Youngblood & Hunter (1999). This line is significant only between 38 dex and 39.9 dex and 
fits our data points over this range too. The slope of the our long dashed line and that 
of Youngblood & Hunter (1999) are similar. The diamonds above the four leftmost points 
are the incompleteness corrected values. Note that they do not change the shape of the 
histogram much. 



Fig. 13. — Ha image with the positions of the H II regions in our sample superimposed 
as white squares. These H II regions were chosen by Hllphot to have a signal-to-noise of 
100 or better. Roughly 1000 H II regions were detected at this signal-to-noise threshold. 
Three annuli are drawn on this image to represent where our subsamples are located. Each 
subsample has the same number of H II regions. These subsamples have the same luminosity 
function and size distribution slopes with respect to one another, indicating no change by 
proximity to the super star clusters. The surface density of H II regions is 1.33 times higher 
in the first annuli (the one nearest the super star clusters) than that from the next annuli 
out. 
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Table I. Observational Parameters of HST/WFPC2 New & Archival Data 



Filter 


Band/ 


PI 


GO-Program 


Date 


Exposure Time 




Emission Line 




Number 




(sec) 


F487N 


R(3 


Shopbell 


8133 


23 Sept. 1999 


3200 


F502N 


[0 m] 


Shopbell 


8133 


23 Sept. 1999 


1500 


F547M 


Stromgren y 


Shopbell 


8133 


23 Sept. 1999 


60 


F656N 


Ha 


Shopbell 


8133 


23 Sept. 1999 


1600 


F671N 


[S II] 


Shopbell 


8133 


23 Sept. 1999 


3000 
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Table 2. Measured and Derived Quantities of NGC 1569 (3x3 5cr results) 



Property Units Value 



Fhq 


ergs s 1 cm 2 


3.4xl0~ n 




ergs s _1 


1.9xl0 40 


SFRhq- 


M yr- 1 


0.15 


Qo(Ha) 


photons s _1 


1.4xl0 52 


FHa,non— pi 


ergs s -1 cm -2 


5.0xl0~ 12 


LHa,non— pi 


ergs s _1 


2.9xl0 39 


LHa,non-pi/LHa,total 




0.15 


Anon— pi/Atotal 




0.16 


L mechanical 


ergs s _1 


0.56-1.1x10 


LfJa, mechanical /^-'Ha, total 




0.1-0.2 



Note. - 



"non-pi" indicates non-photoionized results. 
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Table 3. H II Size Distribution and Luminosity Function Equations 



Property Units 

Diameter Lower & Upper Limits pc 
Size Distribution Equation 
Adjusted R 2 

Luminosity Lower & Upper Limits ergs s" 
Luminosity Function Equation 
Adjusted R 2 

Luminosity Lower & Upper Limits ergs s" 
Luminosity Function Equation 
Adjusted R 2 



Value 



6.3-44.8 

log(N)=-3.02±0.27xlog(D[pc])+5.16±0.27 
-0.9 

^q36.81_^q38.91 

log(N)=-l±0.08xlog(L[ergs s- 1 ])+39.25±2.14 
-0.92 

■^g37.22_^g38.91 

log(N)=-1.22±0.06xlog(L[ergs s- 1 ])+47.5±2.14 

-0.97 



Note. — la uncertainties come after the ±. Adjusted R 2 values are a measure of how well the 
linear regression fits the data given their uncertainties. The higher the number (maxing at 1), the 
better that the fit is. Significant regressions usually have values between 0.9 and I. A negative 
value indicates an inverse correlation. 
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